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SUMMARY 


The  DNA  Underground  Technology  Program  (UTP)  is  developing  a  method  to  predict 
structural  hardness  for  use  in  survivability  and  vulnerability  assessments  of  underground 
facilities.  This  method  is  proceeding  through  a  combination  of  theoretical,  analytical,  and 
experimental  activities.  As  part  of  the  UTP  experimental  program,  DNA  contracted 
RE/SPEC  Inc.  to  perform  laboratory-scale  tests  on  intact  sedimentary  rock.  This  work  is 
motivated  by  the  need  to  analyze  underground  structures  probabilistically  to  account  for  the 
large  amount  of  scatter  in  mechanical  properties  exhibited  by  rock  and  the  uncertainty  m 
site  characterization  and  imprecision  in  loads.  The  overall  objectives  of  this  laboratory 

tj»gtir>g  program  are  to 

1.  Develop  a  constitutive  model  for  intact  sedimentary  rock  (i.e.,  Salem  Limestone)  that 
is  valid  for  the  large  strain,  large  shear  stress,  low  pressure  (<150  MPa  confining 
pressure),  and  complex  strain  histories  expected  in  the  rock  immediately  surrounding 
deep  tunnels  subjected  to  ground-shock  loading;  and  to, 

2.  Construct  a  database  adequate  for  DNA  to  develop  a  probabilistic  constitutive  model 
for  intact  sedimentary  rock. 

Marginal  probability  density  functions  and  a  correlation  matrix  were  determined  for  the 
parameters  of  the  Weidlinger  cap  plasticity  model.  This  information  is  sufficient  to 
transform  non-normal  correlated  random  variables  into  statistically  independent  normal 
ones.  The  experimental  program  comprised  precision  laboratory  testing,  parameter 
estimation,  and  probability  distribution  model  fitting.  A  novel  testing  scheme  allowed  a 
complete  set  of  material  parameters  to  be  determined  from  a  single  test  specimen.  Statisti¬ 
cal  properties  and  distributions  were  determined  by  replication.  The  fitted  model  was 
validated  by  predicting  loading  histories  not  included  in  the  fitting  database.  The  validation 
results  clearly  demonstrated  that  the  parameter  estimation  procedure  produced  high-quality 
parameter  values,  and  that  the  fitted  model  could  not  only  reproduce  measured  data 
included  in  the  fitting  database,  but  also  that  it  could  predict  responses  from  tests  conducted 
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SECTION  1 


INTRODUCTION 

ffigh-quality  solutions  to  complex  underground  structural  mechanics  problems  require 
that  constitutive  relations  be  sufficiently  realistic  and  that  uncertainties  be  understood  and 
quantified.  There  are  several  types  of  uncertainties.  Uncertainty  can  be  categorized  as  (1) 
inherent  (variability  inherent  in  the  material  or  environment),  (2)  statistical  (incompleteness 
of  statistical  data),  (3)  modeling  (use  of  simplified  models),  and  (4)  human  (errors,  quahty, 
judgment).  Only  inherent  uncertainty  is  irreducible.  A  probabilistic  approach  can  be  used  to 
treat  both  irreducible  and  reducible  uncertainties.  Irreducible  uncertainties  can  be  modeled 
as  random  variables.  Reducible  uncertainties  can  be  treated  by  assuming  that  the  mean  and 
standard  deviation  of  the  reducible  random  variables  are  random,  which  results  in  a 

bound  to  the  calculated  probability.  In  this  report  we  deal  only  with  inherent 
material  variability  to  account  for  the  large  amount  of  scatter  in  the  mechanical  properties 
exhibited  by  rock.  Constitutive  models  for  porous  sedimentary  materials  are  reasonably  well 
developed  and  can  simulate  a  variety  of  phenomena  such  as  extension  fractinre  in  tension, 
shear  failure  under  deviatoric  stresses,  and  pore  collapse  under  compressive  mean  stresses. 
Cap  plasticity  models  represent  this  type  of  behavior.  Constitutive  equations  under  the 
categoiy  of  cap  plasticity  models  have  been  formulated  by  numerous  authors,  e.g.,  Drucker 
et  al.  (1957),  DiMaggio  and  Sandler  (1971),  Miller  and  Cheatham  (1972),  Baron  et  al.  (1973), 
Baladi  and  Hadala  (1974),  Sandler  et  al.  (1976),  Sandler  (1976),  Sandler  and  Rubin  (1979), 
Yamada  and  Abou-Sayed  (1979),  Baladi  and  Rohani  (1979),  Celle  and  Cheatham  (1981), 
Cheatham  et  al.  (1984),  Katona  (1984),  Katona  and  Mulert  (1984),  Resende  and  Martin 
(1985),  McCarron  and  Chen  (1987),  Chu  and  Brandt  (1987),  Simo  et  al.  (1988),  Huang  and 
Chen  (1990),  Rubin  (1990),  Abduljauwad  et  al.  (1992),  Hofstetter  et  al.  (1993),  and  Schwer 
and  Murray  (1994).  Cap  models  developed  by  these  authors  describe  material  behavior 
ranging  from  very  simple  to  veiy  complex.  Some  of  the  more  complex  features  include 
material  behavior  at  finite  strain,  anisotropy,  rate  dependence,  nonlinear  isotropic  and 
kinematic  hardening  of  shear  and  volumetric  initial  yield  surfaces  (compacting  and  dilating), 
pore  pressure,  and  Lode-angle  dependence  of  yield  surfaces. 


Procedures  for  cap  model  fitting  are  described  by  several  of  the  authors  cited.  Laborato¬ 
ry  tests  used  most  often  include  uniaxial  strain,  triaxial  compression,  and  hydrostatic 
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compression  tests.  Some  of  the  more  sophisticated  material  models  require  cyclic  stress- 
strain  and  high  strain-rate  tests. 

In  the  traditional  approach,  material  parameters  are  evaluated  sequentially.  For 
example,  the  unloading  portion  of  a  hydrostatic  compression  test  (slope  of  pressure  versus 
volumetric-strain  curve)  is  used  to  determine  the  bulk  modulus;  the  Tinloading  portion  of  a 
triaxial  compression  test  (slope  of  stress-difference  versus  strain-difference  curve)  is  used  to 
determine  the  shear  modulus;  a  series  of  proportional  loading  or  constant  mean  stress 
triaxial  compression  tests  is  used  to  estabUsh  a  shear  yield  envelope.  Parameter  estimation 
procedures  fit  a  shear  yield  (or  failure)  function  to  these  data.  With  elastic  and  shear  yield 
parameters  held  fixed,  the  remaining  parameters  of  the  cap  portion  of  the  model  are 
determined  by  an  iterative  (trial-and-error)  simulation  procedure;  i.e.,  a  cap  shape  and 
hardening  rule  are  postulated  and  the  behavior  of  the  model  calculated  and  compared  to 
measured  material  data  from  uniaxial  strain,  triaxial  compression,  or  hydrostatic  compres¬ 
sion  tests.  This  type  of  fitting  approach  gives  approximate  physical  interpretation  to  some  or 
all  of  the  parameters. 

Nevertheless,  at  least  three  deficiencies  plague  the  traditional  approach.  First,  the  data 
matrix  requires  a  relatively  high  number  of  laboratory  tests  and  therefore  a  high  number  of 
test  specimens.  For  example,  a  three-parameter  shear  yield  function  requires  at  least  three 
triaxial  compression  or  proportional  loading  tests  to  determine  the  parameters  uniquely. 
However,  many  more  than  three  tests  are  usually  needed  because  of  data  scatter.  The 
Americsm  Society  for  Testing  and  Materials  (ASTM  Standard  D  2664, 1989)  recommends 
that  at  least  nine  tests  be  conducted  at  different  confining  pressures  to  determine  a  single 
set  of  shear-strength  properties.  Similar  arguments  hold  for  tests  required  to  evaluate  the 
remaining  parameters.  A  second  weakness  is  that  errors  in  values  of  the  first  parameters 
evaluated  propagate  into  values  of  subsequent  parameters.  Parameter  values  may  therefore 
depend  on  the  order  in  which  they  are  determined.  The  third  disadvantage  is  that  the 
traditional  approach  does  not  permit  evaluation  of  the  statistical  correlation  that  may  exist 
among  parameters. 

These  issues  pose  significant  problems  when  seeking  material-parameter  probability 
density  functions.  If  the  traditional  approach  were  used  as  the  estimation  procedure  to 
compute  point  estimates  of  parameters  and  replication  used  to  determine  statistical 
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properties,  the  whole  series  of  experiments  would  have  to  be  repeated  many  times  and  the 
estimation  procedure  applied  to  each  data  set  in  turn  to  form  a  large  sample  drawn  from  the 
sampling  distribution.  Such  a  procedure  would  clearly  include  a  prohibitive  number  of  tests, 
requiring  hundreds  of  specimens. 

The  report  summarizes  recent  efforts  to  determine  marginal  probability  distributions 
and  the  correlation  among  material  parameters  of  a  Salem  limestone  cap  plasticity  model. 
The  cap  model  is  the  one  proposed  originally  by  DiMaggio  and  Sandler  (1971),  generalized  by 
Sandler  et  al.  (1976),  and  modified  by  Sandler  and  Rubin  (1979).  The  authors  devised  a 
t^ating  scheme  and  parameter  estimation  procedure  allowing  all  the  material  parameters  to 
be  determined  from  a  single  specimen.  Comparisons  of  predictions  against  measurements 
from  other  test  types  (e.g.  uniaxial  strain  tests)  and  from  tests  performed  under  loading 
histories  different  from  those  found  in  the  fitting  database  (e.g.  triaxial  compression  tests 
conducted  at  different  confining  pressures)  validate  the  procedure.  Probability  density 
functions  are  determined  from  twenty  replicate  tests  providing  data  to  construct  histograms 
for  each  cap  model  parameter.  The  Kolmogorov-Smimov  (K-S)  test  (see  Ang  and  Tang, 

1975)  measures  relative  goodness*of-fit  among  candidate  probability  distribution  fimctions. 
Distributions  are  identified  that  best  characterize  the  parameter  histograms.  Also  deter¬ 
mined  is  the  matrix  of  correlation  coefficients  giving  the  degree  of  correlation  among  the 
parameters. 

From  a  practical  standpoint,  results  of  this  study  will  be  key  in  underground  probabil¬ 
istic  structural  analyses.  While  it  is  generally  conceded  that  a  full  Monte  Carlo  solution  is 
impractical  for  large-scale  underground  calculations,  considerable  progress  has  been  made  in 
the  development  and  use  of  approximate  analytical  methods.  Notable  among  these  are 
reliability  analysis  methods  for  geotechnical  applications  based  on  Most  Probable  Point 
Locus  (MPPL)  concepts  in  which  Fast  Probability  Integration  (FPI)  is  used  to  estimate  the 
most  probable  combination  of  random  variables  that  lead  to  failure  (for  examples  in  the  use 
of  this  technique  for  geotechnical  applications  see  Harren  and  Fossum,  1991;  Fossum  et  al., 
1991;  and  Fossum  and  Munson,  1994).  Advanced  mean  value  iteration  techniques  compute 
converged  MPPL  and  estimate  probability  of  failmre.  These  methods  are  acciuate  relative  to 
Monte  Carlo  results  even  in  the  tails  of  the  distributions.  Where  a  Monte  Carlo  analysis 
might  require  100,000  deterministic  solutions,  the  FPI  method  might  require  as  few  as  50 
deterministic  solutions.  It  is  important  to  note,  however,  that  these  FPI  methods  require 
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that  the  random  variables  be  normally  distributed  and  statistically  independent.  If  the 
random  variables,  such  as  the  constitutive  parameters,  are  non-normal  correlated  variables, 
they  must  be  transformed  to  statistically  independent  normal  variables.  This  can  be 
achieved  if  marginal  distributions  and  the  correlation  matrix  for  the  parameters  are  known. 
The  results  obtained  in  this  report  are  sufficient  for  this  purpose. 

The  remainder  of  the  report  describes  the  cap  plasticity  model,  experimental  program, 
and  the  procedures  used  to  evaluate  the  parameters,  their  statistical  correlations,  and 
distributions.  A  discussion  of  resffits  follows  including  those  of  a  validation  study. 
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SECTION  2 


CAP  PLASTICITY  MODEL 

A  detailed  description  and  a  solution  algorithm  for  the  cap  plasticity  model  can  be  found 
in  Sandler  and  Rubin  (1979).  Because  the  imposed  test  conditions  were  known  precisely,  it 
was  possible  to  construct  a  response  model  corresponding  to  each  test.  All  of  the  tests  were 
conducted  in  an  axially-symmetric  compression  machine  giving  rise  to  triaxial  compression 
stress  states.  Using  a  negative  sign  convention  for  compression  and  denoting  the  maximum, 
intermediate,  and  minimum  principal  stresses  as  Oi,  Oj,  O3  respectively,  a  state  of  triaxial 

compression  is  said  to  exist  when 

Oi  =  O2  >  O3 

A  lateral  pressure  in  the  testing  machine  corresponds  to  Oi  =  O2 ,  while  the 

imposed  axial  stress  is  equal  to  O3.  The  square  root  of  the  second  invariant  of  deviator 
stress,  in  this  system,  becomes 

^  ^  <2.2) 

where  the  stress  difference,  Ao,  is  the  difference  between  the  confining  pressure  and  the 
imposed  axial  stress,  Ao  =  -  O3.  The  first  stress  invariant,  A,  is  defined  by 

A  =  o,  +  02  4  03  =  3P  <2-3) 


where  P  denotes  “mean  pressure.” 

The  cap  model  (Sandler  and  Rubin,  1979)  was  formulated  in  terms  of  two  invariants,  I, 
and  Ja-  A  three  invariant  model  is  generally  required  for  rocks  and  simUar  brittle  materials 
such  as  concrete.  If  a  model  with  more  than  two  invariants  is  used,  these  invariants  can 
also  be  expressed  in  terms  of  the  imposed  test  conditions.  However,  depending  upon  the 
functional  form  of  such  a  model,  additional  tests  may  be  needed  (but  not  necessarily)  in  more 
than  just  triaxial  compression  stress  states  to  evaluate  all  of  the  material  parameters.  In 
our  study,  we  used  Equations  2.2  and  2.3  to  develop  response  models  in  terms  of  the 
imposed  test  conditions  since  aU  data  were  obtained  at  the  same  Lode  angle,  which 
corresponds  to  triaxial  compression. 
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The  model  can  be  used  to  model  elastic  behavior,  shear  failure,  and  cap  plasticity.  The 
model  is  a  classical  rate-independent,  associative  plasticity  model  in  which  the  yield  surface 
consists  of  a  shear  yield  surface  and  a  hardening  cap  as  illustrated  in  Figure  2-1.  In  our 
work  we  have  chosen  to  fix  the  shear  yield  surface.  While  this  surface  has  been  called  a 
shear  failure  surface,  it  is  more  correctly  termed  a  limit-state  or  ultimate-shear-yield 
surface.  Holding  this  surface  fixed  is  not  a  necessary  model  restriction.  If  experimental 
data  warrant,  a  hardening  shear  surface  can  be  included  as  part  of  the  complete  model. 
Linear  elastic  behavior  occurs  when  the  stress  is  within  the  composite  shear-failure  and 
yield-cap  surfaces.  The  two  elastic  constants  determined  firom  the  test  data  are  the  bulk 
modulus,  K,  and  the  shear  modulus,  G. 

When  the  stress  point  lies  on  the  shear-failure  envelope,  shear  failure  occurs  according 
to  the  criterion 

^  =  A  -  CexpiSBP)  (2.4) 

V3 

where  A,  B,  and  C  are  material  constants.  During  shear  failure,  the  plastic  strain  comprises 
a  gbpar  component  and  a  volumetric  (dilatant)  component. 

When  the  stress  point  lies  on  the  cap  and  pushes  it  outward,  plastic  strain  comprises  an 
irreversible  decrease  in  volume  called  compaction,  and  a  shear  component  for  other  than 
pure  hydrostatic  compression  stress  states.  The  cap  motion  is  related  to  the  plastic  decrease 
in  volume  through  a  hardening  rule.  The  shape  of  the  cap  is  described  as  an  elliptical 
surface  defined  by 

^  -  Uk)?  -  [3J>  -  LMf  (2.6) 

V3  « 

in  which 

i'"’  "  {  0  if  t  2  0  ■  -BW  -  Cexp(S«]  «.6) 

The  cap  position  parameters  LCk)  and  X(k)  locate  the  current  cap  surface.  The  material 
parameter,  R,  defines  the  ratio  of  principal  ellipse  radii  of  the  cap  surface.  The  hardening 
parameter,  k,  is  defined  through  a  functional  of  XCk)  and  volumetric  plastic  strain,  , 
caused  only  by  cap  action. 
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Figure  2-1.  Ultimate  shear  and  cap  yield  surfaces  of  the  cap  plasticity  model. 


in  which  W  and  D  are  material  parameters,  Xq  is  the  initial  cap  position,  and  is  a  history- 
dependent  functional  of  given  by 

Tp“  /\\  (2.8) 

Sandler  and  Rubin  (1979)  introduced  Equation  2.8  to  prevent  the  cap  from  retracting. 

An  important  issue  that  wiU  not  be  addressed  here  because  it  does  not  affect  parameter 
estimation  is  the  manner  in  which  plastic  strain  rate  is  determined  when  loading  develops 
on  both  cap  and  shear  yield  surfaces.  This  occurs  in  the  current  model  at  the  point  of 
horizontal  tangency  to  the  cap  eUipse  in  Ao  -  P  space  where  the  cap  intersects  the  shear 
yield  surface.  For  a  more  detailed  discussion  of  this  issue  the  interested  reader  should 
consult  articles  by  Resende  and  Martin  (1985)  and  Simo  et  al.  (1988).  Load  paths  in  the 
current  work  are  chosen  so  that  simultaneous  loading  does  not  occur  on  both  yield  surfaces. 


SECTION  3 


EXPERIMENTAL  PROGRAM 


3.1  SPECIMENS. 

Salem  limestone,  also  called  Spergen,  Indiana,  and  Bedford  limestone,  is  a  porous, 
bioclastic,  brecciated,  calcite-cemented,  medium-strength  limestone  (ASTM,  1992).  Its 
carbonate  phases  consist  of  69  percent  fossiliferous  calcite  and  31  percent  calcite  cement. 
Porosity  varies  between  about  12  percent  and  18  percent  depending  on  quarry  location.  The 
block  sample  used  for  specimen  preparation  had  a  density  of  2.698  Mg/m®  and  a  porosity  of 
13.2  percent. 

Specimens  of  this  rock  are  fabricated  in  the  form  of  solid  right-circular  cylinders  having 
nominal  diameters  of  55  mm  and  lengths  of  138  mm  for  length-to-diameter  ratios,  L.D,  of 
2.5.  The  specimens  are  first  cored  from  large  blocks  using  a  conventional  rock  coring  barrel 
moimted  in  a  vertical  milling  machine.  Then  the  ends  of  the  cores  are  finished  to  final 
ijimpnginna  by  Cutting  each  specimen  to  approximate  length  in  a  slabbing  saw  and  trimming 
the  ends  in  a  lathe.  Both  the  coring  and  cutting  operations  use  water  as  a  coolant. 
Following  the  two-step  preparation  procediire,  each  specimen  is  subjected  to  an  acceptance 
test.  The  acceptance  test  includes  specifications  for  flatness  (<  0.02  mm),  straightness  (<  0.3 
mm  gapping),  and  perpendicularity  (<  0.001  radians).  These  specifications  are  more 
stringent  than  those  given  by  the  ASTM  Standard  D4543-85  (ASTM,  1988)  and  were 
selected  to  isolate  only  material  variability  by  reducing  specimen  dimension  uncertainties 
and  piiminating  eccentric  loading  problems  for  poorly  machined  specimens.  In  addition  to 
these  stringent  dimensional  tolerances,  the  maximum  allowable  difference  between  any 
measured  diameter  and  the  average  of  six  measurements  taken  for  three  sets  of  orthogonal 
diameters  at  the  mid-height  and  near  the  specimen  ends  is  not  allowed  to  exceed  0.1  mm. 
Likewise,  the  length  is  not  allowed  to  vary  by  more  than  1.0  mm  firom  the  required  L:D 
multiplied  by  the  average  diameter.  Specimens  failing  the  acceptance  test  are  discarded  or 
remachined  until  each  specification  is  met. 

Specimens  are  soaked  in  ordinary  tap  water  for  30  minutes  to  raise  their  water 
saturation  level  to  50  percent.  Each  specimen  is  then  placed  between  two  stainless  steel 
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metal  endcaps  and  the  endcap-specimen  interfaces  covered  with  lengths  of  1.75-mm-thick 
Viton  tubing  that  extend  to  within  6  mm  of  the  mid-height  of  the  specimen.  Two  thin,  0.3- 
mm-thick  jackets  are  placed  over  the  entire  specimen-endcap  assembly  and  sealed  to  the 
endcaps  with  lockwire.  The  inner  and  outer  jackets  are  fabricated  from  Viton  and  latex, 
respectively.  The  advantages  of  this  jacketing  technique  are  summarized  in  the  next  section. 

3.2  LABORATORY  TEST  SYSTEMS. 

The  specimens  are  tested  using  a  computer-controlled,  servohydraulic  test  system,  Model 
315.03  manufactured  by  MTS  Systems,  Inc.  The  frame  and  axial  load  actuator  are  rated  at 
an  n-gifll  load  capacity  of  -4.5  MN  in  compression.  The  pressure  vessel  and  pressurizing 
systems  used  for  triaxial  compression  testing  are  rated  at  confining  pressimes  of  -150  MPa. 
The  axial  force  anH  confining  pressure  systems  operate  independently  through  the  use  of  two 
separate  servo  loops  that  increase,  decrease,  or  maintain  loads  according  to  a  command 
signal  from  a  computer  that  controls  the  test. 

Because  mechanical  testing  provides  data  to  quantify  imcertainties  in  material  parame¬ 
ters  and  because  these  parameters  are  used  to  quantify  uncertainties  in  structural  behavior, 
it  is  clear  that  test-induced  imcertainties  should  be  small.  The  authors  endeavor  to 
minimize  the  imcertainty  in  the  data  caused  by  the  test  system  through  use  of  (1)  optimal 
transducer  location,  (2)  automated  test  control,  and  (3)  appropriate  calibration  techniques. 

Transducer  location  is  very  important  to  these  tests.  The  test  machine  is  equipped  with 
a  load  cell  inside  the  pressure  vessel  and  extensometry  that  is  moimted  directly  to  the 
specimen.  Frictional  forces  are  eliminated  from  the  load  cell  reading  by  placing  the  load  cell 
inside  the  vessel  as  shown  in  Figure  3-la.  Additionally,  the  cell  is  compensated  for  confining 
pressure  so  that  its  output  is  nominally  zero  at  hydrostatic  stress  states.  Even  at  high 
hydrostatic  pressures,  a  high  gain  can  be  applied  to  the  load  cell  signal  to  improve  resolution 
of  the  load  measmement  and  enhance  precise  axial  load  control. 

Specimen  deformation  is  measmed  using  direct-contact  extensometers.  Axial  and 
circumferential  direct-contact  extensometers  are  mounted  to  the  specimen  over  the  jackets 
(Figure  3-lb).  The  circiimferential  extensometer  is  placed  between  the  ends  of  a  roller-link 
f'hain  wrapped  about  the  mid-height  of  the  specimen  and  measures  changes  in  chord  length 
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defined  by  the  ends  of  the  chain.  The  initial  chord  length  is  typically  about  20  mm.  The 
chain  is  positioned  over  the  two  thin  jackets  but  does  not  contact  the  thick  Viton  tubing  that 
covers  the  specimen-endcap  interfaces.  Mellegard  et  al.  (1993)  have  shown  that  jacket  type 
and  thickness  have  a  profound  effect  on  deformation  measurements  of  the  circumferential 
extensometer.  Therefore,  extensometer  measurements  are  corrected  for  jacket  deformations 
based  on  calibrations  performed  on  6061-T6  aluminum.  The  axial  extensometer  has  a  50.8- 
mm  gage  length  which  is  set  to  within  0.0125  mm  using  pins  in  the  opposing  arms  of  the 
extensometer.  Each  set  of  opposing  arms  is  positioned  on  either  side  of  the  specimen  so  that 
the  gage  length  of  the  extensometer  is  centered  over  the  midheight  of  the  specimen.  Contact 
between  the  extensometer  and  specimen  occurs  through  pins  located  in  the  arms  of  the 
extensometer  and  held  in  place  with  springs.  Measurements  by  the  axial  extensometer  are 
unaffected  by  deformation  of  the  protective  jacket.  The  benefit  in  using  extensometers  that 
mount  directly  to  the  specimen  as  shown  in  Figure  3-lb  is  that  deformations  are  measured 
within  the  zone  of  rock  that  is  presumed  to  have  a  uniform  strain  field.  Further,  any  non¬ 
specimen  deformations  of  the  platens,  specimen-platen  interfaces,  or  machine  distortion  are 
not  part  of  the  deformation  measurement.  As  with  the  load  cell,  mounting  the  extensometer 
directly  on  the  specimen  enhances  test  control  because  of  the  increased  extensometer  signal 
sensitivity  to  specimen  deformation. 

Use  of  servo-hydraulic  testing  techniques,  implemented  and  controlled  by  a  prepro¬ 
grammed  computer,  greatly  reduces  imcertainty  in  the  data.  One  of  the  main  advantages  in 
using  a  computer  is  repeatability  of  loading  conditions  between  specimens.  With  a  non¬ 
computerized  system,  repeatability  of  loading  is  limited  by  operator  skill.  With  a  computer¬ 
ized  system,  the  only  limit  is  resolution  of  the  digital  device  used  to  produce  loading 
commands.  The  system  used  here  employs  a  12-bit  digital-to-analog  converter  to  produce 
commands  that  drive  the  servo-control  loops.  This  provides  control  resolution  of  1  part  in 
4,096,  which  means  that  for  a  full-scale  confining  pressure  range  of  -150  MPa,  the  pressure- 
level  command  has  a  precision  of  approximately  0.037  MPa. 

Of  the  three  actions  taken  to  reduce  imcertainty  caused  by  the  test  system,  the  most 
conspicuous  one  involves  calibration.  Calibration  in  this  work  follows  the  ANSI/ASME 
Performance  Test  Code  on  Measurement  Uncertainty  (PTC  19.1)  that  states  that  insofar  as 
possible  “the  calibration  process  should  include  a  reasonable  simulation  of  instrument  test- 
Uke  conditions.”  Thus,  during  calibration,  transducers  are  connected  in  their  normal 
orientations  on  the  test  system  and  their  outputs  are  recorded  through  the  analog-to-digital 
converter  at  the  computer.  A  typical  calibration  consists  of  applying  twenty  known  standard 
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(a)  Load  frame 


Figure  3-1.  Schematic  of  load  frame  and  deformation  transducers  used  for  all  testing. 
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inputs  to  the  transducer  and  reading  the  corresponding  transducer  outputs  at  the  data 
collection  point.  The  correlation  between  the  transducer  outputs  and  known  standard  inputs 
provides  the  sensitivity  and  offset  for  that  transducer.  Standard  inputs  are  provided  by 
standards  that  are  traceable  to  the  U.S.  National  Institute  for  Standards  and  Technology 
(NIST).  The  analog-to-digital  converter  uses  13  bits  to  digitize  the  magnitude  of  the  signal 
and  one  bit  to  indicate  sign.  Thus,  measurements  are  made  with  a  precision  of  1  part  in 
8,192  over  full  scale  griH  are  not  limited  by  resolution.  Measurement  accuracies,  however, 
are  limited  by  the  standards  used  for  calibration  and  repeatabihty  of  the  transducers. 
Calibration  verification  is  performed  by  determining  the  difference  between  indicated 
readings  versus  standard  input  after  applying  the  standard  input  in  10  equal  steps  over  the 
calibrated  range.  The  transducer  response  at  each  step  is  predicted  using  previously 
determined  cahbration  constants.  The  process  of  cahbration/verification  is  acceptable  if 
predicted  values  fall  within  1.0  percent  of  reading  for  each  step  in  the  verification. 

The  typical  calibration  method  described  above  is  appropriate  for  load  and  confining 
pressure  transducers.  The  cahbrated  ranges  of  these  transducers  are  -2.5  MN  and  -150 
MPa,  respectively.  The  range  on  load  permits  measurements  of  an  absolute  stress  difference 
to  1,100  MPa  for  specimens  with  nominal  diameters  of  55  mm.  The  ranges  in  load  and 
pressure  encompass  the  ranges  estabhshed  by  the  test  matrix.  The  typical  cahbration 
method  is  also  appropriate  for  deformation  transducers,  but  was  modified  as  described  below 
to  improve  resolution  and  accuracy. 

The  stress  ranges  required  by  the  test  matrix  cause  large  strains  that  pose  a  unique 
experimental  problem  when  used  in  model  fitting  efforts  aimed  at  estimating  both  elastic 
and  plastic  model  parameters.  That  problem  is  how  to  maintain  high  resolution  and 
accuracy  for  measurements  of  smaU  elastic  deformations  while  simultaneously  aUowing  for 
the  measurements  of  large  plastic  deformations.  The  problem  is  addressed  by  first  increas¬ 
ing  the  gain  of  the  signal  conditioner  that  services  each  deformation  extensometer  so  that  a 
fuh-scale  voltage  of  10  volts  represents  approximately  one-tenth  of  the  physical  range  of  each 
extensometer.  This  action  produces  an  extensometer  output  signal  that  is  highly  sensitive  to 
smaU  extensometer  measurements  and  very  accurate  over  a  smaU  range  of  specimen 
deformation.  Then,  to  utUize  the  full  physical  ranges  of  the  extensometers  for  accommodat¬ 
ing  large  plastic  strains,  a  unique  feature  of  the  test  system  is  employed  that  apphes  an 
adjustable  offset  voltage  to  each  signal  received  fi*om  the  deformation  extensometers.  This 
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adjustable  offset  voltage  is  controlled  by  the  system  computer  and  is  set  to  shift  the  output 
of  the  signal  conditioner  back  to  zero  whenever  the  level  approaches  10  volts.  The  computer 
keeps  track  of  the  current  offset  voltage  used  to  shift  the  signal  and  applies  appropriate 
calibration  factors  when  converting  the  data  acquisition  signal  into  engineering  imits. 
Multiple  shifts  of  the  siEmal  will  occur  to  accommodate  specimen  deformations  that  approach 
the  full  physical  ranges  of  the  extensometers. 

The  technique  of  amplifying  the  extensometer  signal  and  then  offsetting  the  signal 
conditioner  output  provides  superb  sensitivity  in  test  control  and  greater  resolution  is 
realized  in  data  acquisition  at  all  levels  of  deformation  and  loading.  This  approach,  however, 
requires  a  more  complex  calibration  of  the  extensometers  than  the  typical  calibration 
described  previously.  A  typical  transducer  calibration  involves  the  determmation  of  only  two 
calibration  constants;  i.e.,  the  sensitivity  and  offset.  The  sensitivity  and  offset  represent  the 
slope  and  intercept  of  the  line  fitted  to  calibration  data  pairs  comprising  standard  inputs 
versus  indicated  readings.  The  calibration  procedure  for  the  extensometers,  however, 
requires  the  determination  of  fovir  calibration  constants;  two  that  represent  the  nominal 
slope  and  intercept  as  determined  in  a  typical  calibration,  and  two  more  that  represent  how 
the  nnminfll  slope  changes  whenever  the  signal  is  shifted.  The  calibration  procediue  begins 
by  setting  the  extensometer  to  its  nat\iral  null  position  that  represents  the  begiiming  of  its 
operating  range.  At  this  point,  the  signal  conditioner  output  will  be  at  zero  without  any 
need  for  external  shifting  of  the  signal  level.  With  the  signal  conditioner  gain  set  for  high 
amplification,  a  typical  calibration  run  is  made  over  the  small  range  of  deformation  allowed 
by  the  high  amplification.  At  the  end  of  the  first  typical  calibration,  the  extensometer  is 
positioned  at  one-tenth  of  its  physical  range  and  the  signal  conditioner  output  is  10  volts. 
This  output  voltage  is  then  shifted  back  to  zero  by  applying  the  external  voltage  and  another 
typical  cahbration  nm  is  now  made  with  the  “null”  extensometer  position  shifted.  These 
steps  of  calibration  and  voltage  shifting  continue  until  the  physical  range  of  the 
extensometer  is  reached.  This  proced;ire  of  calibrating  over  small  strain  ranges  that  are 
concatenated  by  external  voltage  shifts  to  form  a  large  strain  range  provides  enough  data  to 
determine  aU  four  calibration  constants.  The  full  ranges  of  the  axial  and  circximferential 
extensometers  can  handle  large  plastic  strains  to  -10  percent  and  6  percent,  respectively. 

Although  direct-contact  transducer  measurements  are  preferable  to  other  types  of 
measxirements,  their  use  in  the  harsh  environment  of  the  pressure  vessel  of  the  test  system 
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raises  questions  about  the  effect  of  pressure  and  temperature  on  the  output  of  transducers. 
Senseny  (1987)  has  shown  that  pressure  affects  only  the  offset  of  the  extensometers  while 
temperature  affects  both  offset  and  sensitivity.  The  effects  of  pressure  on  the  deformation 
transducers  used  in  this  study  were  examined  and  reported  by  Mellegard  et  al.  (1993). 

Using  their  results,  the  changes  in  the  offsets  of  the  axial  and  circumferential  extensometers 
are  -0.15x10'*  m/MPa  and  0.42x10'*  m/MPa,  respectively.  The  change  in  the  offset  of  the 
internal  load  cell  because  of  pressure  is  -9.3x10'®  MN/MPa.  Although  all  testing  is  performed 
at  a  constant  temperature  of  20*C  in  this  study,  temperature  changes  of  the  confining 
pressure  fluid  and  therefore  the  direct-contact  transducers  can  occur  as  a  result  of  adiabatic 
heating  of  the  ^nnfiniTigr  fluid  during  pressurization.  Mellegard  et  al.  (1993)  also  studied  this 
phenomenon  to  examine  the  effect  of  pressurization  rates  on  the  outputs  of  direct-contact 
transducers  and  concluded  that  the  effects  are  small  provided  absolute  pressurization  rates 
are  less  tiian  or  equal  to  0.02  MPa/s.  Therefore,  in  this  study,  hydrostatic  loading  stages  are 
performed  at  a  pressmization  rate  of  -0.02  MPa/s. 

3.3  EXPERIMENTAL  PROGRAM. 

The  load  history  imposed  allows  the  material  parameters  to  be  determined  from  a  single 
test  specimen  with  little  xmcertainty  and  generally  with  minimal  confoimding  among  other 
parameters.  The  load  path  includes  several  stages.  Because  test  conditions  in  each  loading 
stage  are  known  precisely,  response  models  can  be  constructed  corresponding  to  each  stage. 
Nonlinear  regression  can  then  be  used  to  evaluate  the  parameters  such  that  the  fitted  model 
gives  the  best  fit  to  the  data  in  a  least  squares  sense.  Response  models  refer  to  predicted 
strain  in  stress-rate  control  and  to  predicted  stress  in  total  strain-rate  control.  Simulations 
of  tests  follow  exactly  the  load  path  imposed  experimentally. 

The  first  stage  of  the  test  is  hydrostatic  compression  in  stress-rate  control.  In  this  stage, 
hydrostatic  pressure  is  applied  to  a  cylindrical  test  specimen  at  a  constant  rate  of  -0.02 
MPa/s  imtil  -150  MPa  is  reached.  This  value  is  at  or  near  the  value  needed  to  initiate 
plastic  volume  compaction  (pore  collapse).  Axial  and  radial  strains  are  measured  during  the 
test. 

The  second  stage,  triaxial  compression  loading,  begins  with  a  mode  switch  from  stress- 
rate  control  to  strain-rate  control.  Radial  and  tangential  stresses,  identical  from  symmetry. 


16 


are  held  constant  at  -150  MPa  while  the  axial  strain  is  controlled  at  a  constant  rate  of  -1.0  x 
lO-^/s  \mtil  an  axial  strain  of  -8.5  percent  is  reached.  During  this  stage  radial  strain  and 
axial  stress  are  measured.  This  stage  is  terminated  before  the  shear  surface  is  reached  but 
only  after  a  m^or  portion  of  the  cap  surface  has  been  traversed. 

The  third  stage  is  unloading  with  a  mode  switch  back  to  stress-rate  control.  Eadial  and 
tangential  stresses  remain  constant  at  -150  MPa  while  axial  stress  is  reduced  to  -150  MPa  at 
a  constant  rate.  Axial  and  radial  strains  are  measured  during  this  stage.  Figure  3-2  in  Aa- 
P  space  shows  the  load  paths  corresponding  to  these  three  stages.  Figure  3-2  also  shows  the 
position  of  the  yield  surfaces  relative  to  the  load  paths.  After  these  first  three  stages,  the 
measured  stress-strain  data  provide  sufficient  information  to  evaluate  sensitivity  of 
simulated  responses  to  changes  in  elastic  and  cap  parameters. 

To  obtain  data  sufficient  to  define  the  shear  yield  function  parameters,  the  next  portion 
of  the  test  includes  a  series  of  hydrostatic-unload/triaxial-load  stages  that  we  term  a 
“walkdown.”  The  test  specimen  is  unloaded  hydrostatically  firom  the  hydrostatic  stress  state 
that  existed  after  the  third  stage  (-150  MPa).  When  a  hydrostatic  stress  state  of  -40  MPa  is 
reached,  the  specimen  is  loaded  in  axial  strain-rate  control  vmtil  the  axial  stress  stops 
increasing.  The  specimen  is  then  unloaded  to  a  hydrostatic  stress  level  of  -40  MPa.  The 
specimen  is  then  unloaded  hydrostatically  to  -35  MPa  and  the  process  repeated.  Five  data 
points  are  collected  during  this  “walkdown”  procedure  in  -5  MPa  hydrostatic  increments. 

Stress-strain  data  from  the  first  three  stages  and  the  “walkdown”  strength  data  provide 
sufficient  information,  from  a  single  specimen,  to  determine  all  nine  parameters  of  the  cap 
plasticity  model.  Statistical  properties  and  sampling  distributions  for  the  material  parame¬ 
ters  follow  from  replicate  testing  with  20  specimens. 

The  20  rephcate  tests  are  part  of  a  larger  testing  program  conducted  to  characterize  the 
rate-independent  deformation  of  Salem  limestone.  The  larger  program  comprises  87 
hydrostatic/triaxial  compression  tests  and  10  uniaxial  strain  tests.  In  addition  to  the  20 
tests  performed  at  a  confining  pressure  of  -150  MPa,  the  hydrostatic/triaxial  compression 
test  matrix  includes  20  unconfined  tests,  10  tests  each  at  -50  MPa  and  -100  MPa,  and  one 
test  at  each  5  MPa  ^nnfining  pressure  increment  between  -5  MPa  and  -145  MPa.  These 
tests  and  the  \iniaxial  strain  tests  provide  data  for  model  validation  and  are  not  used  in  the 
parameter  estimation  procedure. 
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Figure  3-2.  Three-stage  load  paths  and  their  positions  relative  to  yield  surfaces. 


3.4  TEST  PROCEDURE. 


An  instrumented  specimen  is  placed  inside  the  pressure  vessel  of  the  test  system.  The 
pressure  vessel  is  then  filled  with  silicone  oil  and  a  small  compressive  hydrostatic  preload 
(<1  MPa)  is  manually  applied  to  the  specimen  by  the  system  operator.  Control  of  the  test 
system  is  then  given  over  to  the  computer.  The  control  software  requests  operator  input 
such  as  specimen  dimensions,  load  path  specifications,  etc.,  and  then  begins  to  load  the 
specimen.  For  hydrostatic/triaxial  compression  tests,  the  test  system  initially  applies 
hydrostatic  stress  at  a  rate  of -0.02  MPa/s  (Stage  1,  Figure  3-2).  When  the  target  hydrostat¬ 
ic  stress  is  reached,  the  axial-load  servo  loop  switches  to  axial  strain  control  and  axial 
loading  continues  at  constant  confining  pressure  at  an  axial  strain  rate  of  -IxlO-^/s  (Stage  2, 
Figure  3-2).  When  approximately  85  percent  of  the  physical  range  of  the  axial  extensometer 
has  been  traversed,  the  axial-load  servo  loop  switches  back  to  load  control  and  the  stress 
difference  and  the  hydrostatic  stress  are  removed  sequentially  (Stage  3,  Figure  3-2).  For 
.miavigl  strain  tests,  axial  load  is  again  applied  in  axial  strain  control;  however,  lateral 
straining  is  constrained  by  increasing  confining  pressure.  The  axial  extensometer  is  used  as 
feedback  for  the  axial  load  servo  loop,  while  the  circumferential  extensometer  is  used  as 
feedback  for  the  confining  presstire  servo  loop.  The  imposed  axial  strain  rate  is  -1x10  /s, 
which  is  one  order  of  magnitude  lower  than  used  in  the  hydrostatic/triaxial  compression 
tests.  This  lower  strain  rate  is  used  to  prevent  high  confining  pressxire  rates  that  would 
result  in  adiabatic  heating  of  the  direct-contact  extensometers.  Uniaxial  strain  loadmg 
continues  until  the  fuU  range  of  the  confining  pressure  transducer  is  traversed  at  which  time 
the  algebraic  sign  of  the  axial  strain  rate  is  reversed  and  the  specimen  is  unloaded. 

During  both  the  hydrostatic/triaxial  compression  tests  and  the  umaxial  strain  tests,  data 
are  acquired  based  on  the  output  levels  of  both  the  axial  extensometer  and  the  internal  load 
fp]]  When  the  output  of  either  of  the  transducers  changes  by  more  than  0.01  percent  of  full 
scale  (based  on  physical  range),  data  are  logged.  Logged  data  include  time,  axial  load, 
confining  pressure,  axial  deformation,  and  circumferential  deformation.  A  typical  datafile  for 
a  hydrostatic/triaxial  compression  test  contains  about  1,000  data  points. 


Data  reduction  is  performed  using  a  commercial  spreadsheet  software  package.  In  this 
study,  straius  are  true,  or  logarithmic  strains,  and  stresses  are  true,  or  Cauchy  stresses;  i.e.. 
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the  instantaneous  dimensions  of  the  specimen  as  determined  from  the  deformation  measure¬ 
ments  are  used  to  calculate  stress. 

In  the  model  fitting  procedure,  the  number  of  data  points  for  each  test  is  reduced  from 
about  1,000  points  to  30  points  for  each  of  the  three  loading  stages  (hydrostatic  compression, 
triaxial  compression  loading,  and  triaxial  compression  unloading)  for  a  total  of  90  points  per 
specimen.  The  number  of  points  is  reduced  for  reasons  of  practicality  and  to  place  equal 
weight  on  each  stage  of  the  load  path.  Data  points  are  selected  at  equally  spaced  intervals 
within  the  original  datafile  using  the  algorithm,  I  =  (N  ■  l)/29  where  I  is  the  interval 
between  selected  data  points  and  N  is  the  total  number  of  data  points  acqiiired  for  a  given 
stage  of  the  load  path.  Therefore,  if  1,000  points  are  contained  within  the  original  datafile, 
the  interval  between  selected  data  points  is  34.  The  first  and  last  data  point  of  each  stage 
are  always  included  in  the  reduced  datafile. 

3.5  EXPERIMENTAL  RESULTS. 

Figures  3-3  and  3-4  show  the  stress-strain  data  from  all  20  tests.  Axial  stress  is  plotted 
against  total  axial,  lateral,  and  volumetric  strains.  These  figiu-es  illustrate  the  relatively 
small  variability  in  the  material.  The  inelastic  strain  in  these  tests  is  clearly  the  result  of 
voliune  compaction  as  shown  in  the  unloading  portion  of  the  stress-volumetric  strain  cmves. 
Figure  3-5  shows  peak  values  of  Ag/^  plotted  against  3P  for  the  “walkdown”  testing  stages 
of  the  tests.  The  data  indicate  that  the  shear  limit  smface  is  concave  toward  the  3P  axis. 
Figure  3-6,  using  data  from  one  of  the  test  specimens,  illustrates  the  type  of  stress-strain 
behavior  exhibited  by  each  test  specimen  dmdng  the  “walkdoAvn”  stages  of  the  tests.  From 
these  stress-strain  plots  one  can  clearly  see  the  relatively  ductile  nature  of  this  limestone. 


20 


Axial  Stress  (MPa) 


Strain 


Figure  3-3.  Axial  stress  versus  axial  strain  and  radial  strain  for  aU  hydrostatic/triaxial 
compression  tests. 
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Axial  Stress  (MPa) 


Volumetric  Strain 


Figure  3-4.  Axial  stress  versus  voliunetric  strain  for  all  hydrostatic/triaxial  compression  tests. 
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Stress  Difference  (MPa) 
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Figure  3-6.  Stress-strain  behavior  during  the  “walkdown”  stages  of  the  ultimate  shear  tests. 


SECTION  4 


PARAMETER  ESTIMATION/STATISTICAL  PROPERTIES/DISTRIBUTIONS 


4.1  MODEL  FITTING  PROCEDURE. 


In  the  global  nonlinear  regression  analysis,  data  from  all  stages  of  a  single-specimen  test 
are  fitted  simiiltaneously  using  response  models  appropriate  for  the  load  path  followed  in 
each  testing  stage.  Different  response  models  are  required  for  each  stage  because  different 
quantities  are  controlled  during  each  stage.  The  objective  function  minimized  in  fitting  the 
cap  model  is  the  weighted  s\un  of  squared  residuals  (RSS)  given  by 


RSS  = 


i-l  }  STAGE  1 

30  SO 

E^l/H^ssD®  (ok  -  0*33/  + 

j-1  i-l 

30  30 

(ek  -  ek)*  E(l/»«33l)^  (ek  ' 

••1  i-l 

6E(1/I1Ao|1)^  (Aak  - 


STAGE  2 


STAGE  3 


(4.1) 


1-1 


SHEAR  STAGES 


in  which  the  subscript  m  denotes  “measured”  and  the  subscript  p  denotes  “predicted.” 

Terms  enclosed  in  double  vertical  lines  are  maximum  measured  values  for  the  variable 
indicated  and  serve  as  weighting  functions. 

Total  volumetric  strain,  is  the  dependent  variable  in  Stage  1;  axial  stress,  033,  and 
radial  strain,  e^,  are  the  dependent  variables  in  Stage  2;  radial  strain,  Eh,  and  axial  strain, 
£33,  are  the  dependent  variables  in  Stage  3;  and  Aa  is  the  dependent  variable  in  the 
ultimate  shear  stages.  A  weighted  least  squares  estimation  rather  than  a  standard  least 
squares  estimation  is  used  because  error  variance  is  not  homogeneous.  Weighting  also 
accounts  for  different  units  in  the  response  model  and  thereby  normalizes  all  test  informa- 
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tion.  In  addition,  since  there  are  only  five  data  points  in  each  shear  stage,  the  shear  stages 
are  multiplied  by  6  so  that  each  loading  stage  is  weighted  equally. 

The  Bio  Medical  Data  Processing  (BMDP)  statistical  software  package,  (BMDP  Statisti¬ 
cal  Software,  Inc.,  1992),  is  used  to  solve  for  the  response  and  to  fit  the  model  to  laboratory 
data.  A  5th-order  Runge-Kutta  integration  method  is  used  together  with  a  pseudo-Gauss- 
Newton  algorithm  to  perform  nonlinear  least  squares  fitting.  The  regression  analysis  gives 
estimates  of  each  parameter,  standard  error  for  each  parameter,  and  the  correlation  matrix 
of  parameter  estimates. 

4.2  PARAMETER  ESTIMATES. 

The  weighted  least  squares  routine  is  applied  to  the  objective  fimction  defined  by 
Equation  4.1  for  each  of  the  20  single  specimen  tests.  Parameter  estimates  for  each 
specimen  are  given  in  Table  4-1.  The  average  coefficient  of  variation  for  the  parameters  is 
30  percent.  The  highest  correlation  between  parameters  exists  between  A  and  C  with  a 
correlation  coefficient  near  1.0,  The  elastic  parameters  show  the  least  correlation  with  other 
parameters.  In  general,  a  correlation  coefficient  near  ±  1.0  between  two  different  parame¬ 
ters  suggests  that  the  set  of  fitted  parameters  may  be  unstable  and  that  a  condition  of 
multicoUinearity  may  exist  so  that  these  parameter  values  may  not  be  imique.  Two  causes 
of  multicoUinearity  are  the  foUowing:  (1)  the  model  may  inherently  include  parameters  that 
depend  on  one  another,  and  (2)  the  database  may  not  be  sufficient  to  evaluate  certain  model 
parameters  uniquely.  To  reduce  multicoUinearity,  one  could  (c)  modiiy  the  functional  form 
of  the  model,  or  (6)  expand  the  database  of  measxirements.  Nonetheless,  because  the 
correlation  matrices  are  invertible,  it  is  possible  to  transform  the  non-normal  correlated 
parameters  into  statisticaUy  independent  ones,  (e.g.,  Liu  and  Der  Kiimeghian,  1986)  for  use 
in  probabiUstic  analyses. 

Figures  4-1  and  4-2  show,  respectively,  the  model-fits  to  single-specimen  stress-strain 
data  with  the  lowest  and  highest  RSS  values,  while  Figimes  4-3  and  4-4  show  the  corre¬ 
sponding  model-fits  to  ultimate  shear  yield  data  from  the  same  specimens.  The  best  model- 
fit  has  an  RSS  value  of  0.0403,  whUe  the  worst  model-fit  has  an  RSS  value  of  0.0979.  The 
median  RSS  value  for  the  twenty  specimens  is  0.0637. 
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Table  4-1.  Parameter 


Specimen 

Number 

K(MPa) 

G(MPa) 

ACMPa) 

B(l/MPa) 

E-03 

SL1214 

14,562 

8,613 

209.89 

1.810 

SL1231 

18,472 

10,878 

217.04 

1.623 

SL1232 

16,710 

9,593 

220.29 

1.471 

SL1241 

16,432 

8,941 

209.04 

1.648 

SL1242 

15,679 

8,893 

212.86 

1.696 

SL1243 

16,824 

9,766 

186.28 

1.868 

SL1244 

15,802 

10,168 

205.08 

1.683 

SL1245 

15,181 

10,143 

193.58 

1.869 

SL1246 

14,401 

8,558 

186.29 

2.018 

SL1251 

16,805 

9,882 

210.23 

1.801 

SL1253 

17,902 

9,854 

196.15 

1.777 

SL1255 

16,768 

9,259 

209.61 

1.787 

SL1261 

17,378 

10,802 

216.36 

1.712 

SL1265 

16,846 

9,399 

213.01 

1.723 

SL1266 

17,357 

9,461 

197.26 

2.107 

SL1275 

13,671 

8,670 

154.38 

2.998 

SL1276 

16,071 

9,094 

191.12 

2.120 

SL1278 

13,627 

8,414 

185.44 

2.244 

SL1713 

17,385 

9,674 

203.67 

1.947 

SL1716 

16,611 

10,559 

171.06 

2.648 

for  each  specimen. 


C(MPa) 


198.80 


190.07 


195.76 


190.52 


201.69 


166.45 


186. 


175.83 


171.83 


198. 


181.33 


198.49 


201.72 


196.30 


189.54 


152.92 


183.41 


181.96 


190.08 


163.84 


D(l/MPa) 

E-04 


R 

X,(MPa) 

RSS 

E-02 

5.88 

-397.27 

7.41 

6.18 

-442.80 

6.67 

5.90 

-428.16 

4.99 

6.00 

-427.69 

7.08 

5.98 

-422.01 

5.63 

6.67 

-426.55 

7.42 

5.98 

-443.13 

8.97 

5.03 

-601.87 

8.25 

5.19 

-441.51 

6.41 

5.91 

-427.43 

5.92 

6.35 

-431.88 

8.35 

5.63 

-442.56 

9.79 

5.96 

-430.80 

6.77 

5.03 

-480.29 

4.03 

5.49 

-441.87 

5.49 

5.34 

-425.37 

6.72 

5.10 

-427.81 

4.79 

6.27 

-439.47 

8.99 

6.06 

-482.03 

5.30 

5.66 

-464.82 

6.32 

Axial  Stress  (MPa) 


Figure  4-1.  Lowest  RSS  model-fit  to  single-specimen  stress-strain  data. 
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Figure  4-2.  Highest  RSS  model-fit  to  single-specimen  stress-strain  data. 
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Figure  4-3.  Lowest  RSS  model-fit  to  ultimate  shear  yield  data  from  a  single  specimen. 
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Figure  4-4.  Highest  RSS  model-fit  to  ultimate  shear  yield  data  from  a  single  specimen. 


4.3  EVALUATION  OF  STATISTICAL  PROPERTIES  AND  DISTRIBUTIONS. 


To  facilitate  selection  of  probability  density  distribution  functions,  parameter  estimates 
from  the  20  tests  are  plotted  as  a  frequency  distribution  in  a  graphic  representation  known 
as  a  histogram.  The  histogram  breaks  the  total  range  into  a  series  of  intervals  and  indicates 
the  numbers  of  occurrences  into  each  interval  by  rectangles  of  different  sizes.  Representa¬ 
tive  of  such  a  plot  is  the  one  for  parameter  W  shown  in  Figure  4-5.  In  this  case,  indications 
are  that  functions  skewed  toward  higher  values  of  the  random  variable  are  appropriate. 
Distributions  capable  of  capturing  the  skewed  nature  of  the  parameter  are  selected  as 
candidates  including  lognormal,  gamma,  and  WeibuU  (Type  HI  -  extreme  value  distribution). 
For  comparison,  we  also  include  normal  and  uniform  symmetric  distributions.  To  discrimi¬ 
nate  among  candidate  functions,  a  statistic  is  needed  to  evaluate  relative  goodness-of-fit  of 
each  function  against  observed  data.  Two  common  statistics  used  for  this  p\u*pose  are  the 
Chi-squared  (x^)  and  Kolmogorov-Smimov  (K-S)  tests  for  goodness-of-fit  (Ang  and  Tang, 
1975).  The  x*  statistic  is  a  test  on  differences  between  hypothesized  and  observed  frequen¬ 
cies  within  each  interval  of  the  histogram.  The  degrees  of  fireedom  in  a  x^  test  depend  on  the 
number  of  intervals  used  to  create  the  histogram.  This  means  that  the  significance  level  for 
testing  sixitability  of  a  proposed  distribution  function  can  vary  depending  on  how  the  interval 
is  chosen  (Benjamin  and  Cornell,  1970).  In  addition,  it  is  commonly  recommended  that  at 
least  five  observations  be  recorded  within  each  interval.  The  x*  test  results  become  critically 
dependent  on  these  factors  for  small  sample  sizes  and  conflicting  results  can  occur  even 
when  the  same  data  set  is  used  to  test  a  hypothesis  and  only  the  choice  of  histogram  interval 
is  varied. 

An  alternative  to  the  x^  test  is  the  K-S  test  that  compares  differences  between  hypothe¬ 
sized  grirf  observed  cumulative  distribution  functions  (CDF’s).  The  degrees  of  freedom  in  the 
K-S  test  is  unique  and  depends  only  on  the  number  of  data  points  used  in  the  test.  Because 
the  number  of  parameter  estimates  is  relatively  small  (20),  the  K-S  test  is  chosen  as  the 
statistic  to  determine  the  relative  goodness-of-fit  of  the  candidate  probability  distribution 
functions.  The  largest  difference  between  the  experimental  and  predicted  CDF’s  represents 
the  K-S  statistic,  D. 

The  candidate  probability  density  distribution  functions  are  fitted  to  the  histogram  data 
depicted  in  Figure  4-5.  Table  4-2  gives  the  K-S  goodness-of-fit  statistic  for  each  of  the 
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Density 
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Figure  4-5.  Histogram  and  probability  function  for  parameter  W. 
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Table  4-2.  Candidate  probability  density  functions  and  D 
statistic  for  parameter  W. 


CDF 

D 

Normal 

0.255 

Lognormal 

0.218 

Gamma 

0.239 

Weibull 

0.233 

Uniform 

0.421 

candidate  probability  density  functions  for  parameter  W.  The  table  shows  that  the  cmnula- 
tive  distribution  function  corresponding  to  the  lognormal  probability  density  function  gives 
the  lowest  value  of  D.  An  illustration  of  how  well  the  lognormal  CDF  fits  the  experimental 
CDF  for  W  is  shown  in  Figure  4-6.  Figure  4-5  shows  the  lognormal  probability  density 
function  superimposed  on  the  histogram  of  parameter  W.  Distributions  for  the  remaining 
parameters  are  determiued  in  a  similar  fashion.  A  summary  of  these  results  is  presented  in 
Table  4-3  along  with  values  for  the  mean,  median,  and  standard  deviation  for  each  parame¬ 
ter. 
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Comparison  between  predicted  and  experimental  CDF  for  parameter  W. 


Table  4-3.  Parameter  distributions,  statistics,  and  best-fit  values. 


Parameter 

Probability 

Density 

Fimction 

Mean 

Value 

Standard 

Deviation 

Median 

Value 

Best-Fit 

Value 

K(MPa) 

WeibuU 

16,214 

1,358 

16,610 

16,046 

G(MPa) 

Lognormal 

9,531 

740 

9,527 

10,093 

A  (MPa) 

Weibijll 

199.43 

16.80 

204.32 

269.52 

B  (1/MPa) 

Lognormal 

1.922E-03 

3.680E-04 

1.806E-03 

1.211E-03 

C(MPa) 

WeibuU 

185.77 

13.69 

189.81 

251.19 

D  (1/MPa) 

WeibiiU 

8.616E-04 

5.055E-04 

7.437E-04 

6.410E-04 

W 

Lognormal 

0.108 

0.064 

0.074 

.109 

R 

WeibiiU 

5.63 

0.41 

5.65 

5.35 

Xo(MPa) 

Lognormal 

-441.26 

24.20 

-435.67 

-458.60 
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SECTION  5 


BEST  FIT  AND  MODEL  VALIDATION 

5.1  OVERALL  BEST-FIT  PARAMETER  ESTIMATES. 

The  best  estimators  of  the  true  values  of  a  set  of  parameters  are  those  for  which  the 
fitted  model  gives  the  best  fit,  in  a  least  squares  sense,  to  the  data  considered  simultaneous¬ 
ly  fi-om  all  tests.  That  is,  the  objective  fimction  minimized  in  fitting  the  cap  model  includes 
twenty  contributions,  each  one  with  a  form  identical  to  Equation  4.1.  This  set  of  parame¬ 
ters,  shown  in  the  last  column  of  Table  4-3,  is  the  best  set  to  use  in  a  deterministic  analysis 
and  are  commonly  called  the  “best-fit”  parameters.  From  the  overall  fit,  parameter  Xq  has 
the  lowest  mean  coefficient  of  variation  at  5.48  percent,  while  parameter  W  has  the  highest 
value  at  59.26  percent.  Statistically,  this  means  that  parameters  and  W  have  the  least 
and  most  uncertainty,  respectively.  Use  is  made  of  this  set  of  parameters  in  the  validation 

study. 

The  correlation  matrix  presented  in  Table  5-1  is  determined  from  the  overall  best-fit 
calculations  and  is  the  one  to  use  in  probabilistic  analyses  that  make  use  of  the  marginal 
probability  distributions  for  the  cap  model  parameters  determined  in  Section  4. 

5.2  VALIDATION  AGAINST  TESTS  OUTSIDE  THE  FITTING  DATABASE. 

Because  of  the  unconventional  nature  of  the  parameter  estimation  procedure,  a 
validation  study  is  used  to  assess  the  predictive  capability  of  the  fitted  model.  The  overall 
best-fit  parameter  values  are  used  to  predict  responses  for  tests  and  test  t3q)es  not  included 
in  the  fitting  database.  These  tests  include  triaxial  compression  tests  at  different  constant 
confining  pressures,  uniaxial  strain  tests,  and  independent  ultimate-shear  tests  conducted 
with  multiple  test  specimens. 

The  first  validation  problem  predicts  hydrostatic/triaxial  compression  tests  in  which  the 
specimen  is  loaded  to  -100  MPa  in  hydrostatic  compression  followed  by  axial  straining  to  -8.5 
percent  in  strain  rate  control.  The  confining  pressure  remains  constant  at  -100  MPa.  This 
test  is  identical  to  the  first  two  stages  of  the  tests  discussed  earlier  except  that  the  confining 


37 


Table  5-1.  Matrix  of  correlation  coefficients  giving  the  degree  of  correlation 
among  parameters. 


K 

G 

A 

B 

C 

D 

W 

R 

Xo 

K 

1.0000 

G 

-0.0284 

1.000 

— 

0.0894 

-0.0419 

1.000 

(S5Tnmetric) 

B 

-0.2610 

0.0053 

-0.9537 

1.000 

C 

0.0706 

-0.0442 

0.9997 

-0.9511 

1.000 

D 

-0.0869 

-0.2806 

-0.6133 

0.6406 

-0.6166 

1.000 

W 

0.1161 

0.2967 

0.5996 

-0.6270 

0.6019 

-0.9980 

1.000 

R 

-0.0679 

-0.0785 

-0.3943 

0.2266 

-0.3980 

0.3028 

-0.3001 

1.000 

Xo 

0.2618 

-0.1251 

0.7181 

-0.8686 

0.7157 

-0.5423 

0.5230 

0.2270 

1.000 

pressure  at  the  end  of  the  first  stage  is  -100  MPa  rather  than  -150  MPa.  Figure  5-1  shows 
measured  versus  predicted  stress-strain  behavior.  Figures  5-2  and  5-3  show  measured 
versxis  predicted  results  for  similar  tests  in  which  the  first  stage  ends  at  -50  MPa  and  0  MPa 
respectively. 

The  next  validation  problem  predicts  the  results  of  imiaxial  strain  tests.  In  this  test  the 
confining  pressure  is  increased  so  that  the  radial  strain  is  constrained  to  be  zero  while  the 
avifll  strain  increases.  Figure  5-4  shows  the  results  of  ten  uniaxial  strain  tests  versus  the 
predicted  stress-strain  behavior  made  with  the  best-fit  parameters  determined  previously. 
Note  that  the  prediction  appears  to  give  a  best-fit  to  the  uniaxial  strain  data  even  though 
this  test  type  was  not  included  in  the  fitting  database. 

The  filial  validation  exercise  involves  ultimate  shear  yield  tests  in  which  a  different 
specimen  is  used  for  each  of  the  different  confining  pressxires  used  to  define  the  yield 
envelope.  The  measured  versus  predicted  results  are  shown  in  Figure  5-5.  Each  of  the 
measured  data  points  is  taken  from  a  different  specimen.  The  predicted  resiilts  make  use  of 
the  best-fit  parameters  listed  in  Table  4-3  determined  from  replicated  tests  that  included  the 
“walk-down”  load  histories. 
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Axial  Stress  (MPa) 


Predicted  versus  measured  stress-strain  behavior  for  the  hydrostatic/triaxial  compression  validation  tests 
conducted  at  a  confining  pressure  of  -100  MPa. 


Axial  Stress  (MPa) 


Figure  5-2.  Predicted  versus  measured  stress-strain  behavior  for  the  hydrostatic/triaxial  compression  validation  tests 
conducted  at  a  confining  pressure  of  -50  MPa. 


Axial  Stress  (MPa) 
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conducted  at  zero  confining  pressure. 


Axial  Stress  (MPa) 


Figure  5-4.  Predicted  versus  measured  stress-strain  behavior  for  the  uniaxial  strain  validation  tests. 


Aa/ V3(MPa) 
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Figure  5-5.  Predicted  versus  measured  values  of  Ao/^S  plotted  against  3P  for  the  independent  ultimate  shear  tests. 


The  validation  results  clearly  demonstrate  that  the  parameter  estimation  procedure 
produces  high-quality  parameter  values,  and  that  the  fitted  model  can  not  only  reproduce 
measured  data  included  in  the  fitting  database,  but  also  predict  responses  fi-om  tests 
conducted  under  loading  histories  different  from  those  found  in  the  fitting  database. 
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SECTION  6 

DISCUSSION 

Exceptions!  care  was  taken  throughout  this  investigation  to  quantify  inherent  variability 
in  Salem  limestone  by  minimizing  estimation  error,  hiunan  error,  and  modeling  error. 
Carefully  controlled  laboratory  tests  generate  data  for  parameter  evaluation.  Response 
models  are  constructed  corresponding  to  each  test  and  best  estimates  of  the  model  parame¬ 
ters  are  provided  by  the  best  fit  to  the  data  in  a  least  squares  sense.  No  o  priori  assump¬ 
tions  are  required  regarding  the  material  behavior  in  the  laboratory  tests.  All  parameters 
are  evaluated  simultaneously  to  produce  a  consistent  parameter  set.  By  using  a  single 
specimen  to  determine  a  complete  set  of  cap  model  parameters,  statistical  properties  are 
determined  by  straightforward  replication,  thereby  avoiding  the  complicating  assumptions  of 
simulation  methods  involving  repeated  use  of  data.  Such  methods  may  be  unavoidable, 
however,  as  the  level  of  model  sophistication  increases  (e.g.,  accoimting  for  the  effects  of 
unloading-reloading  using  a  combination  of  isotropic  and  kinematic  plastic  hardening  rules, 
or  simulating  the  change  in  shape  of  the  failure  surface  with  the  mean  stress).  In  such 
cases,  Monte  Carlo  computer  simulation  may  be  needed  to  simulate  a  series  of  experiments. 
This  series  could  then  be  replicated  as  many  times  as  necessary  to  obtain  the  relevant 
properties  of  the  sampling  distribution. 
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SECTION  7 


CONCLUSIONS 

Marginal  probability  distributions  and  a  correlation  matrix  were  presented  for  parame¬ 
ters  of  a  cap  plasticity  material  model.  Point  estimates  of  the  parameters  were  determined 
by  a  testing  scheme  and  parameter  estimation  procedme  in  which  the  parameters  were 
determined  from  a  single  specimen.  Statistical  properties  were  determined  from  replicate 
testing,  and  probability  distributions  were  determined  from  candidate  distributions  that  best 
characterized  parameter  histograms. 
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